Haldane-like Chern insulators in buckled lattices with in-plane magnetic fields 
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Chern insulators, first proposed by Haldane over 20 years ago, has not been observed experi- 
mentally. While there is some expectation that topological insulators with broken time reversal 
symmetry may realise the phase, the mechanism - spin orbit interaction together with magnetic 
exchange - has a completely different origin to Haldane's staggered magnetic field. I show that 
buckled lattices with in-plane magnetic fields are a promising route to realising Haldane-like mod- 
els. As a concrete example, I consider silicene, a honeycomb lattice with out-of-plane sublattice 
anisotropy, in an in-plane magnetic field, and show that it is a Chern insulator insulator, even at 
negligibly small magnetic fields, which is analogous to Haldane's original model, though bears some 
interesting differences as well. 

PACS numbers: 73.43.-f,73.43.Cd, 73.43.Jn 



F.D.M. Haldane's 1988 paper, 'Model for a Quantum 
Hall Effect without Landau Levels: Condensed-Matter 
Realization of the "Parity Anomaly" ' [l[ , was a paradigm 
shifting work which is undoubtedly a precursor to the 
field of topological insulators In fact, topological 

insulators in an external field, or in the presence of ferro- 
magnetic exchange coupling, are thought to be Chern, or 
quantum anomalous Hall (Q AH) insulators Q , of which 
Haldane's model is the seminal example. An ultracold 
gas implementation of Haldane's model was recently pro- 
posed, which utilised laser induced pseudomagnctic fields 
to produce a staggered flux Q . 

Both Haldane's original model, and the more recent 
topological insulators, are based on Dirac cone physics. 
The most straightforward statement of this is in terms 
of linear Dirac equations: two adjoining systems with 
massive (gapped) Dirac cones, with Hamiltonians H± = 
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(k x , ky, ±m) • tr, supports gapless edge states at their in- 
terface due to the change in sign of the mass between the 
two systems. In real systems, this amounts to having an 
effective mass term which is oppositely signed to that of 
the vacuum. At the interface between the system con- 
taining the cone, and the vacuum then, there is a point 
at which m — ¥ 0, and the band gap closes. The quan- 
tum anomalous Hall phase is that where there is a single 
(Mod(2)) zero energy species on the edge of a system. 
The quantum spin Hall phase is that where there are 
two (Mod(2)) counter-propagating edge states, as there 
must be if time reversal symmetry is not broken. 

In this Letter, I show that Haldane-like models can be 
constructed, using the same mechanism, of a staggered 
magnetic flux, as in Haldane's original proposal, by em- 
ploying an in-plane magnetic field incident on a buckled 
two dimensional system. Using a concrete example, I 
construct a Hamiltonian which realises the same chiral 
states as in the original model, and yet does not have the 
same chiral structure of flux accumulation. The orienta- 
tion of the incident field can also induce phase transitions 
between the QAH and ordinary insulating states. 




FIG. 1: Hexagonal unit cells with zero net magnetic flux. 
The upper two unit cells are a buckled honeycomb lattice 
(e.g. silicene) in an in-plane magnetic field. Due to the buck- 
led lattice structure, closed loops of nearest and next-nearest 
neighbor bonds have either a positive or negative flux pass- 
ing through them. In the lower pane is shown Haldane's unit 
cell [l[, where the flux (from a spatially inhomogeneous per- 
pendicular magnetic field) through the outer three loops are 
compensated by an opposing flux through the centre. In each 
case, the direction of positive flux accumulation is indicated 
by the arrows along the bonds. Note that I have only shown 
the bonds along a single sublattice for clarity. 



Silicene is a two dimensional honeycomb lattice, where 
the two sites in the unit cell are offset from each other, in 
the out-of-plane direction. I define the plane parallel to 
the lattice, z = 0, such that one triangular sublattice is at 
z = I, and the other at z = —I. The silicene Hamiltonian 
in an in-plane magnetic field and an out-of-plane electric 
field then, is d: 
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ff = *E c ^+ E 4- c ta+ A I>«cK (i) 

where i is the nearest neighbor hopping integral, tL = 

if exp(i§ f? A ■ dl) is the Peierl's substitution modified 
next nearest neighbor hopping, A = IE Z is the electric 
field multiplied by the out-of-plane buckling / = 0.23A 
The out-of-plane field is not necessary to observe the 
QAH effect, but allows us to make contact with Hal- 
dane's work, and when including spin effects, becomes 
important. I have chosen a gauge such that the flux is 
accounted for on next nearest neighbor bonds. We are 
free to choose any consistent gauge such that the total 
phase accumulated over closed loops encircling regions of 
some flux is always commensurate with that flux. Due 
to the definition that z = defines the lattice plane such 
that the two sublattices are at z = ±Z, the particular 
gauge choice in which nearest neighbor flux accumulation 
vanishes, is the Landau gauge, A = Bz(sm9 7 — cos 8, 0), 
where 9 is the angle of the field relative to the horizon- 
tal direction in Fig.pQ. The Hamiltonian, in momentum 
space, becomes 

H % = ct(d (fc)a + d(k) ■ a)c % , (2) 

where ct = (cr tr J), and in which dn(k) = 

cos((j) / (j>o) J^i 2t' cos(fc-i^), where (j> Q = h/2e, and <p/4>a is 
the integrated gauge field over the bond using the Lan- 
dau gauge, which gives the total flux through the loop 
with one next nearest neighbor and two nearest neigh- 
bor bonds ((f> — Bla, with a the lattice constant), as 
indicated in Fig.pQ. 5[ arc three out of six of the next 
nearest neighbor vectors defined by S' n = R(2irn/3)(1, 0), 
where R(0) rotates the vector by 8, and 
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A - sim>/0 o )2f' sin(fc •£■)), 

where Si are the three nearest neighbor vectors defined 
by 8 n = R(2-im/3)(0, 

This particular choice of gauge makes concrete con- 
tact with Haldane's original choice, as well as being in- 
tuitively and computationally convenient. The system, 
together with the inclusion of flux along the bonds, and 
the comparison with Haldane's model, are all shown in 
Fig- [I]. The chirality, noted by Haldane in his model, 
is the clockwise accumulation of positive flux over next 
nearest neighbor hops. The clockwise accumulation of 
positive flux is clear in the lower panel in Fig.pQ. This 
chirality is absent in the buckled lattices (i.e. the upper 
two unit cells shown in Fig.pQ), where hopping over the 
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FIG. 2: Phase diagram for spinless silicene in an in-plane 
magnetic field (such that 4>/<j)o = i"/2), where A = El. Here 
9 is the angle of the in-plane field. The angles tt/6 and n/2 
correspond to the upper left orientation of the field in Fig.pQ, 
whereas the angle n/3 corresponds to the upper right. 

equivalent loop leads to zero phase accumulation. How- 
ever, a triangular loop comprising one next nearest neigh- 
bor and two nearest neighbor hops acquires a non-zero 
positive phase. For the upper left panel, there are two 
loops which acquire a negative phase on clockwise hop- 
ping, and one which acquires a negative phase. There is, 
therefore, a net chirality. In the upper right however, the 
lower triangular plaquette envelopes no magnetic flux, 
and thus there are just two oppositely oriented chiral 
loops, and thus no net chirality. 

The Hamiltonian Eq. [2] describes a graphene-like sys- 
tem with two low energy Dirac cones at the two inequiv- 
alent K points, K± = (±4tt/3,0). At A,B = 0, the 
system is exactly the graphene system, where the low 
energy Dirac cones are massless. A and B then intro- 
duce separate mass terms. The A mass term is due to 
the on-site energy imbalance between the two sublattices 
due to the electric field Stark effect. The magnitude of 
this gap is 2EI. The second, which I call the 'flux gap', is 
the sublattice anisotropy due to the orbital effect of the 
magnetic field, and its magnitude is 

A = v/3*'|2sin(0/2^ o ) - sin(</>/0 o )| sin(30). (4) 

This gap goes as (^>/</>o) 3 for small fields, and so may be 
tiny for realistic laboratory fields. Later, I will discuss 
ways in which this number can be improved by several 
orders of magnitude in the lab. 

Since the flux gap is ex sin(fc) (i.e. the z-component 
of Eq. pi]), it has opposite sign in the two valleys (the 
two Dirac, or K, points). Therefore, A and A^ add, or 
compete, in the two valleys, as is clear from Fig. [3], such 
that the gap in the two valleys, are Ak ± = A ± A^. 
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FIG. 3: Phase diagram including both spins in silicene. For 
small spin orbit coupling and small electric field, we obtain 
two co-propagating chiral bands. This is equivalent to a v — 2 
quantum Hall effect, or equivalently, a doubled Chern insu- 
lator. For an electric field gap larger than the flux gap, the 
system is a trivial insulator, as also noted in Fig. [2] . For small 
electric field gap, but comparatively large spin orbit coupling, 
silicene is a Kane-Mele topological insulator However, 
when Aso +m > A^,, but \Xso —m\< A$, silicene is a Chern 
insulator with a single chiral band, equivalent to a v = 1 
quantum Hall effect. In principle, the diagonal C = 1 section 
of the phase diagram continues indefinitely, so long as the en- 
ergy scales m ~ Aso <C t, such that the K points are not 
mixed by these terms. 



The quantum anomalous Hall (QAH) phase is char- 
acterised by a Chern number of 1 (Mod(2)), where the 
Chern number is the integral of the Berry curvature over 
the Brillouin zone ■ For our Hamiltonian Eq. pQ, it 
takes the particularly simple form 0, [l2[ 
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dd dd k± sgn(A±) 



dk v 



(5) 



where the second equality is the result for the Hamilto- 
nian Eq. pQ expanded to first order in k near the K± 
points. From the second equality, we suspect, and it is 
confirmed by taking the integral Eq. [3] over the entire 
Brillouin zone, that the Chern number is 



C n = 



A > A 
A A < A. 



(6) 



In Fig. [2] is shown the Chern number phase diagram 
of Hamiltonian Eq. [Tj, as a function of in-plane field 
orientation 8, where 8 = corresponds to a field along 
the zig-zag direction, as shown in the upper right portion 
of Fig. pQ. At A = 0, the Chern number is ±1 for almost 
the entire spectrum of field orientations. This can be 



understood by considering Fig.pQ. For all 8 except irn/3, 
there is a positive flux passing through cither one or two 
plaquettes in the unit cell, and a negative flux passing 
through two or one, such that the total flux is always zero. 
The positive fluxes instigate cyclotron motion around the 
triangular plaquettes in one direction, while the negative 
flux instigates cyclotron motion in the opposite direction. 
Therefore we expect to see two co-propagating modes at 
the edge of the material, and one counter-propagating 
mode, leading to a net chiral current. 

It is interesting to note from Fig. [2], that a topological 
phase transition can be induced by varying the orienta- 
tion of the in-planc magnetic field. It is also worth noting 
that at A = 0, and at zero temperature, the system is a 
quantum anomalous Hall insulator for any magnetic field 
strength. 

The primary result of the current work is Fig. [2], rele- 
vant for buckled lattices. Specifically, that chiral edge 
states can be created, akin to Haldane's model intro- 
duced in 1988, for buckled lattices in an in-plane mag- 
netic field. I now turn my attention to the inclusion of 
spin, and some notes on possible experimental realiza- 
tions of this novel result. 

So far I have completely neglected spin. The Chern 
numbers reported in the phase diagram Fig. [2] for a 
fcrmionic system are per spin. For spin degenerate sys- 
tems such as that considered here, each spin species will 
co-propagate. The Zeeman splitting will not affect the 
Chern number, but if it is larger than the flux gap it will 
move the bulk bands of one spin species across the Fermi 
energy. Experimentally, this is an added complication, 
and will be addressed briefly toward the end of the pa- 
per. Therefore, rather than having a total Chern number 
in Hamiltonian Eq. PQ of 1, for a spinful fcrmionic sys- 
tem, there is an extra factor of 2 for the spin degeneracy, 
and the Chern number is in fact 2, corresponding to a 
v = 2 quantum Hall effect, with 2 filled Landau levels. 

In order to produce a spinful fcrmionic system with a 
total Chern number of 1, spin orbit coupling is required. 
Here we add a spin index to Eq. pQ, and introduce the 
Kane-Mele_spin orbit coupling term between next nearest 
neighbors 
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(7) 



where = ±1 depending on whether a left or right turn 
was required in order to hop between the two sites via 
nearest neighbor bonds. In silicene, Aso ~ 4 meV 11 



The low energy expansion of such a term acts as a mass 
with different signs in all of the sublattice, valley, and 
spin degrees of freedom. Therefore it modifies the gap 
Ak ± such that [T3| Ak± s s = A ± A& ± sXso, where 
s = ± corresponds to spin. When Xso + A > A^, a 
mass inversion occurs for a single spin species, in a sin- 
gle valley in the Brillouin zone. This in turn leads to 
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FIG. 4: Dispersion of energy eigenstates for ribbons of: Sil- 
icene at Aso = (upper), and A = 0, the Haldane model 
(middle) with A = 0.05t, and silicene with Aso = 0.07i, 
A = 0.05t (lower). While the Haldane model has a single 
chiral mode per edge, silicene without spin orbit coupling 
has three modes per edge, per spin, one of which is counter- 
propagating with respect to the other two (and one of which 
cannot be seen as it is too close to the bulk bands). There- 
fore the Chern number in the two systems (per spin) is the 
same. In order to obtain a Chern number of 1 in a spinful 
system, spin orbit coupling must be introduced. In the lower 
pane, Aso has caused a band inversion for one K-point, and 
one pseudo-spin, such that this band becomes topologically 
trivial, while the other pseudo-spin retains a Chern number 
of 1. The other parameters for the silicene systems: t' = O.lt, 
(j>/(j>a = tt/2, and 8 = tt/2. 



the destruction of a single edge state, and thus a total 
Chern number of 1. From the phase diagram Fig. [2], 
the way to understand this is that one species has mass 
> A — Aso, and thus C„ = 1, whereas the other has 
A^ < A + Aso, and thus C to t = 1- For larger spin orbit 
couplings however, a second mass inversion takes place, 
when | A — Aso I > A^, such that the system becomes 
equivalent to the original quantum spin Hall insulator 
[2|, with counter- propagating edge states. These phases 
are all indicated in the spinful phase diagram, Fig. [3J. 

In Fig. [4], zig-zag edged ribbons of silicene with and 
without spin orbit coupling, together with the Haldane 
model are shown. For spin degenerate silicene (i.e. ne- 
glecting Zeeman splitting), we see that there is a single 
band (per edge) with three propagation directions, or two 
turning points, per band. These three propagation direc- 
tions correspond to the three cyclotron orbits discussed 
earlier in the context of Fig.pQ. For the sake of com- 
parison, the Haldane model is also shown in the middle 
pane, which just has a single propagation direction per 
edge. These two figures emphasise the contrast of the 
mechanisms of the Haldane and silicene models. 

The lower panel in Fig. [J] shows the effect of the spin 
orbit coupling on the edge states. As discussed, the spin 
orbit interaction inverts the mass of a single spin species 
at a single valley, thus producing a single anticrossing in 
comparison with the spin degenerate case. This in turn 
leads to a Chern number of for one edge state, and a 
Chern number of 1 for the other. 

Finally, I estimate the magnetic field needed to pro- 
duce a Haldane-like model in a real system. For any finite 
magnetic field, at zero temperature, silicene is a Chern 
insulator. However, in order to obtain a chiral edge state 
per spin, an appreciable flux gap must be formed by the 
magnetic field in order to obtain edge states not lost to 
thermal fluctuations. For silicene, the next nearest neigh- 
bor bonds are separated by sa 4 A, and the vertical buck- 
ling amount is w 0.46 A. Therefore </>/0 o « 2 x 10 -5 T _1 . 
Unfortunately, for small fields, the flux-gap Eq. [J goes 
as 3 , as mentioned earlier, and so an extremely large 
field (~ 100 T) is required to obtain a gap of only 2/j,K . 
This is clearly unrealistic. 

A more feasible scenario can be obtained by inducing 
next nearest neighbor hopping anisotropy. As alluded to 
earlier, this can be achieved by the inclusion of a plane 
polarized laser. This can be understood as the hopping 
parameters being adjusted, in the same manner as by 
the Peierl's substitution, by a time averaged gauge field. 
More concretely, in the regime hu <C , this modifies 
the hopping matrix elements such that t — > Jo(eaE / 'hu)t 
by Floquet theory [l4, 15 1, where E — Eq cos(wt)-E, and 
a is the lattice constant. As an example, this modifies 
the gap equation Eq. [J] such that, for an electric field 
polarized in the plane parallel to 9 — 0. and for small 4> 
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A = V3t'(t>/(l)o\{J Q (eaE /2huj) - J a (eaE Q /huj))sin(39)\, 

(8) 

where Jo is the modified Bessel function of the second 
kind. For a low frequency, high intensity laser, the flux 
gap will increase linearly, and appreciably, with <f>. It is 
clear then that if the next nearest neighbor hopping can 
be modified such that the Bessel function difference in 
Eq. [5] is optimised, a measurable flux gap may become 
feasible. As a concrete example, for a 2.5 THz laser, 
where haj = lOmcV, at an intensity of 5 x 10 7 Vm" 1 , a 



flux gap of pa 1.3i'</>/0o is achievable, which, for a 35 
T field, is w 0.5 K. 

Another route to realising anisotropic next nearest 
neighbor hopping is via uni-axial strain. In this case, 
the flux gap is modified as A^ oc 8t'(j>, where t' — > t' + St' 
along the strain axis. However, the added complication 
of emergent gauge fields must be taken into account in 
this case flij . 

In both of the above examples, the Zceman splitting 
will likely cause bulk bands for one K point to cross the 
Fermi surface. For a g- factor of 2, the Zeeman split- 
ting at 35 T is ~ 20 K, which in the specific example 
mentioned above, is much larger than the flux gap. I 
emphasise however, that the Zeeman splitting does not 
affect the Chern number in any way. Away from the K 
points, the bulk bands will still be gapped in general, 
and so the gapless edge states could feasibly be imaged 
with angle resolved photoemission spectra. Alternatively, 
the edge states could be measured via a nonlocal trans- 
port experiment, using an H-bar type geometry, which is 
precisely the way in which 2D topological insulator edge 
states have been compellingly verified jlij ]. In order to 
obtain a bulk gap however, the Zeeman splitting must 
be less than the flux gap, and so the plaquette area must 
be maximized and/or the next nearest neighbor hopping 
must be large, or a material used where the g-factor is 
suppressed. 

More exotic systems may be employed to maximize 
the staggered flux through the unit cell. For instance, 
bilayer systems can be envisaged where the inter-layer 
gap becomes much larger than that reported here, due 
to the vastly increased interlayer spacing. Another pos- 
sibility is that the "designer Dirac fermions" engineered 
by precise placement of carbon monoxide molecules to 
produce "molecular graphene" |17| could be engineered 
to produce a buckled structure with large out-of-plane 
displacement, and/or strong next nearest neighbor cou- 
pling and anisotropy. 

In conclusion, I have introduced a new conceptual 
framework whereby Haldanc-like models can be realised. 
The key material property in this proposal is a buck- 
led lattice, such that an in-plane magnetic field produces 
a unit cell with zero net flux, yet regions with positive 



and negative flux which individually support cyclotron 
orbits. In this way, the chirality necessary in implement- 
ing a Haldanc-like model can be realised. Using silicene 
as a specific example, I showed that over a range of field 
angles, silicene in an in-plane field is topologically equiv- 
alent to the Haldane model per spin, with chiral edge 
states arising by similar, yet distinct physical processes, 
and that topological phase transitions can be induced by 
rotating the orientation of the field. The physical real- 
ization of a Haldane-like model in a condensed matter 
system would be a remarkable achievement. 
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